Multiple beam interference in a quadrupolar glass fiber 
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Motivated by the recent observation of periodic filter characteristics of an oval-shaped micro-cavity, 
we study the possible interference of multiple beams in the far field of a laser-illuminated quadrupo- 
lar glass fiber. From numerical ray-tracing simulations of the experimental situation we obtain 
the interference-relevant length-difference spectrum and compare it with data extracted from the 
experimental filter results. Our analysis reveals that different polygonal cavity modes being refrac- 
tively output-coupled in the high-curvature region of the fiber contribute to the observed far-field 
interference. 
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Introduction. Optical fibers and cavities have attracted 
a lot of interest in recent years. On the theoretical side, a 
plethora of phenomena related to the interplay of classi- 
cal and quantum chaos is found [0J^] . On the experimen- 
tal side, fibers are used either as (active) lasing fibers 
in microlasers, or as (passive) optical filters which are 
of great technological interest for planar integrated fil- 
ter applications. Planar dielectric ring and disc cavities 
have been used as micron-sized optical filters mainly with 
evanescent light coupling, working with nearly total in- 
ternal reflection. However, evanescent coupling between 
the cavity curved sidewall and the waveguide flat sidewall 
requires a very precise fabrication with a gap spacing in 
the sub-^m range. Therefore, filter techniques using non- 
evanescent coupling which allows gap sizes larger than 
sub-fim are technologically desirable. 

Interestingly, recent experiments |^] using an oval- 
shaped micro-cavity have shown periodic output filter 
characteristics which are potentially useful in the above 
context. The purpose of this paper is to present an anal- 
ysis of the experimental data and numerical ray-tracing 
simulations which allow a theoretical understanding of 
these experimental findings. Our results display a peri- 
odic filter spectrum in the far-field interference of refrac- 
tively output-coupled modes of the cavity, in agreement 
with the experiment. We will show that the interfer- 
ence of multiple beam parts corresponding to polygonal 
round-trip orbits can lead to the observed periodic filter 
characteristics in a narrow window of the far-field angle. 

First we briefly summarize the experiment of Ref. || : 
A passive (non-lasing) quadrupolar fiber of high- 
refraction glass (n — 1.8) is illuminated by a laser beam, 
Fig. 0a. The tunable laser source with wavelengths in the 
670 nm range produces a Gaussian beam with a width 
(spot size) of 30 fim, the cavity axes are 150 and 180 
/ito. The far-field elastic scattering spectrum is mea- 
sured with a linear array detector. It shows filter reso- 
nances as function of incoming wavelength with a good 
peak to background ratio of about 40, Fig. ||a, but only 
under very specific input and output coupling angles 0i iO 
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FIG. 1. a) Schematic experimental setup [[| for the fil- 
ter experiment. Shown is the cross section of the oval-shaped 
fiber. If the input is coming from a broadband source, then 
the output shows periodic filter characteristics, see also Fig. ^. 
f?i iD are the input and output angles, b) Typical orbits con- 
tributing to the far-field response, see text. 



("magic window"). The periodicity of the spectrum is a 
clear sign of the interference nature of the phenomenon; 
in addition the filter peaks display inhomogeneous broad- 
ening which is an indication of multimodc interference. 
The light wavelength A is much smaller than the size of 
the cavity, therefore quantum effects can be assumed to 
be small. 

Analysis of experimental data. If we accept that the 
periodic output characteristics can be interpreted as in- 
terference of (classical) rays, then a length analysis of the 
contributing geometric paths appears suitable. We start 
by defining the amplitude-weighted length distribution 
T(L) for the inferfering rays, 



T(L)= Ai5(L-Li 

paths i 
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where Ai and Li are amplitude (at the detector) and 
optical length (from source to detector) of each path i 
hitting the detector. The interference pattern J{k) as 
function of the vacuum wavevector k = 2ir/\ is given by 
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In the last step we have introduced the length-difference 
spectrum S(l), given by the self-convolution (correlation 
function) of T(L), S(l) = J dLT(L)T(l + L); for discrete 
paths with lengths Li the quantity S(l) will be non-zero 
for I = Li — Lj The length-difference spectrum S(l) 

will be the main quantity of our analysis, it is related by 
Fourier transformation to the observed interference pat- 
tern J(fc); the information about absolute path lengths 
does of course not enter the interference result. 

The Fourier transform S(l) of a representative set of 
experimental data from Ref. [pi is shown in Fig. 0. We 
see large contributions to S(l) at roughly equally spaced 
I values, it is tempting to identify the spacing with the 
path length of one round trip corresponding to a single 
dominating cavity orbit. However, the intensity for larger 
I is rather small, and the broad peaks suggest that more 
than one orbit contributes to S(l). Similar results are 
obtained with other data sets of the experiment pj . 

At this point it is useful to discuss possible interfer- 
ence scenarios. The simplest possibility is interference 
of rays tunneling out of a single stable orbit that is 
traced over and over again. In this case T(L) contains 
equally spaced peaks with monotonically decreasing in- 
tensity (~ t(l — t) 1 ^ 1 where t is the tunneling rate and 
i the number of round trips); a similar decay is then 
also found in the peaks of the difference spectrum S(l). 
For the particular case of the quadrupolar fiber a can- 
didate stable orbit is the "diamond" , here the tunneling 
rate is rather small. From this one would expect a slow, 
monotonous decay of the peak intensities in S(l) — obvi- 
ously this is not in agreement with the experimental data. 
Furthermore, in the present experimental geometry the 
diamond-like orbit cannot be excited by refractive input 
coupling, but only by tunneling - the resulting intensity is 
much lower than for refractive coupling, and is too small 
to account for the experimental observation. Other orbits 
like rectangular/trapezoidal modes or whispering-gallery 
(WG) orbits Q are either unstable or cannot be excited 
with the present non-evanescent coupling. We are there- 
fore lead to consider the interference of rays from multiple 
orbits as explanation for the observations of Ref. || . 

Ray-tracing simulations. We focus a two-dimensional 
geometry representing the cross-section of the fiber used 
in the experiment, see Fig. [l^,. The geometry is de- 
fined by the mean radius R of the cavity (R = 82 [im 
in Ref. ||) and its excentricity e » 0.1. The quadrupole 
boundary is given by r{4>) = i?(l + ecos2</>) in polar co- 
ordinates, and the lengths of the half axes are R(l ± e). 

The incoming beam is discretized into a sufficiently 
high number of equally spaced parallel rays, for simplic- 
ity we employ a rectangular beam profile. The intensity 
fraction of each ray that penetrates into the quadrupole 
is given by the Frcsncl formula, its angle by SnelPs law. 
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FIG. 2. Representative experimental spectrum of Ref. Q. 
a) raw interference data, J(k), plotted as function of the 
wavevector k, 9i = 60°, o = 56°. b) Modulus |] 
of the Fourier transform of J(k) corresponding to the 
length-difference spectrum S(l). The strong decay of the peak 
intensity for lengths I longer than 3 or 4 round trips is a char- 
acteristic feature. 



Quantum effects are negligible due to A <C R. The dy- 
namics of each ray is then governed by the laws of a "Fres- 
nel billiard" [pPJ], i.e., by straight propagation, specu- 
lar reflection at the quadrupolar shaped boundaries, and 
evolution of the intensity according to Fresnel's law for 
reflection and transmission. In particular, for TE polari- 
sation H the Fresnel laws for the reflected (transmitted) 
electric field strength E r (E t ) for a incoming field Ei are 
given by || 
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where n is the refractive index of the cavity material, 
n = 1 outside the cavity, and \ is the boundary an- 
gle of incidence, Fig. |l|. Strictly, these relations hold for 
a planar interface; modified Fresnel formulae for curved 
interfaces can be derived ||], but in the present case of 
a very large size parameter, R 3> A, the curvature ef- 
fect is negligible. We assume perfectly reflecting walls 
for angles of incidence x bigger than the critical angle, 
Xc = arcsin(l/n), and exclude leakage due to quantum 
tunneling. 

In the simulation, we trace each ray of the incom- 
ing beam numerically to construct the interference pat- 
tern (2|j^]. For angles X < Xc we allow for refractive 
escape of the part of the ray that is determined by 
Eq. (^), but follow further the remaining part inside the 
quadrupole until its intensity falls below a threshold of 
10~ 6 of the initial intensity due to subsequent subcritical 
reflections. For the transmitted part, we determine the 
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FIG. 3. Results of the ray-tracing simulation, a) Inten- 
sity histogram showing the distribution of orbit lengths L vs 
their output angle o . The gray scale indicates the intensity 
(black - maximal). The input angle is fixed at Bi = 60°. b) 
Intensity vs. orbit length L at a specific detector position 
(in the "magic window", see text), obtained by integrating 
the above histogram over a narrow interval of output angles, 
58° < o < 62° . If we neglect external path differences this 
quantity is equivalent to the length distribution T(L). 



far field angle of the leaving ray again by Snell's law. 

In Fig. |l]b we show a couple of typical trajectories that 
are found upon scanning of the incoming beam. Due to 
the finite excentricity and the finite beam width we find 
not only WG orbits but a large number of orbits which es- 
cape (as well as enter) around the points of highest curva- 
ture of the quadrupole, as known from the study of asym- 
metric resonance cavities ^M. In particular, we find rays 
that undergo several polygonal-like round trips (in con- 
trast to typical WG orbits they come closer to the center 
of the quadrupole) before their intensity eventually drops 
below the threshold. This process of intensity loss may 
happen at one single reflection or (more likely) upon a 
couple of subsequent reflections/transmissions. Due to 
the preferred escape points near the highest wall curva- 
ture 0], the distribution of the orbit lengths, measured 
from entering the fiber until escape, is sharply peaked at 
integer multiples of half the length of a round trip. Orbits 
with approximately 0.5, 1.5, 2.5 round trips are shown in 
Fig. [l]b - these (and similar longer) orbits contribute to 
the far-field interference for an output angle as chosen in 
Fig. |l|a, i.e., for output points opposite to the incident 
beam. 

Given the experimental findings, we are especially in- 
terested in rays with escape direction within the "magic 
window" . The filter characteristics was observed in the 
far-field total intensity - in this situation the detector 
is placed at a distance large compared to the cavity ra- 
dius. Therefore only rays within a narrow window of the 
output angle O are detected, whereas the exact escape 
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0i = 60°, B a = 60°, and e 
are as in Fig. |^. 



0.1, for details see text. The axes 



position on the quadrupolar boundary is of no impor- 
tance. In contrast, in a near field measurement (usually 
done with a focussing lens) a narrow interval of output 
positions is sampled, with a rather large range of output 
angles. 

We now turn to the ray-tracing simulation results for 
the far-field interference. The primary output of the sim- 
ulation are the ray trajectories, their escape points and 
output angles, and the output (transmission) intensities 
which result from the Fresnel formulae ([|). Fig. ||a shows 
lengths L and output angles o in the form of a inten- 
sity histogram plot ( "Fresnel- weighted" histogram) for a 
fixed input angle 9i — 60°. Here, L are the geometric 
orbit lengths inside the cavity in units of R; the optical 
lengths are found from L — nLR + L cxt + iphasc, where 
L cxt are external length differences arising from the dif- 
ferent input and output points of the interfering rays, 
and Lphasc is determined by the phase shifts that occur 
upon the reflections. Note that the quantity shown in 
Fig. ||a is the output-angle-resolved version of the length 
spectrum T(L) defined in (|l|). For most output angles, 
the intensities^are rather small, and arise primarily from 
short lengths L. However, in a narrow window which cor- 
responds to escape points in the high-curvature region, 
the total output intensity is larger, and short as well as 
longer orbits carry significant weight. 

If we now integrate the above quantity over a small 
range of far-field angles corresponding to the angle range 
convered by the detector, we obtain the length spectrum 
T(L), shown in Fig. ||b. The peaks are easily found to 
correspond to orbits as shown in Fig. [j]b (and longer 
ones). After converting L into optical lengths L, we ob- 
tain from T(L) the length-difference spectrum SU) by 
self-convolution, and by Fourier transformation Q the 
resulting interference pattern J(k) as function of the 
wavevector k = 27r/A, Fig. ||. 
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Qualitative agreement between experiment and theory 
is most easily checked by comparing the length-difference 
spectra S(l), Figs. ||b and[|b. The main feature, namely a 
number of roughly equally spaced peaks with comparable 
intensity, but very little intensity at larger lengths, is 
nicely reproduced by the ray-tracing data. 

We note that one cannot expect to reach a quantita- 
tive agreement with experiment due to the uncertainty 
in experimental parameters like size and excentricity of 
the fiber and angle of the illuminating laser beam. The 
precise interference pattern depends strongly on these pa- 
rameters, in particular it is extremely sensitive to length 
changes of the order of the light wavelength. For this rea- 
son we have neglected both L cxt and L p hasc when plot- 
ting Fig. ||. These contributions to the optical length are 
of the order of one or several wavelengths, and are cer- 
tainly smaller than the uncertainty in the fiber size. We 
have checked that the interference pattern depends only 
slightly on the numerical discretization procedure used 
for the incoming beam. 

By varying input and output angles, our simulation 
data clearly show that in most situations the rays from 
all parts of the incident beam travel in polygonal cavity 
orbits only for a very short time (up to 2 round trips) 
before they leave the cavity via refraction. Only for a 
narrow range of input angles an appreciable part of the 
beam leads to polygonal orbits with a longer lifetime, 
which then are refractively output-coupled after a larger 
number of round trips into a narrow window of far-field 
angles. The far-field output angle window depends sensi- 
tively on the input angle, so we predict that the "magic 
window" as observed in Ref. will move (in the far-field 
angle) with varying input angle of the beam. In partic- 
ular, for 9i = 70° we found the magic window at output 
angles 10° smaller than for B L = 60°, whereas for 8i = 50° 
the filter effect almost disappears. 

In addition, a number of conditions have to be met to 
observe the cavity filter effect: (i) A finite excentricity 
of the fiber is needed to produce (chaotic) orbits which 
come close to the center of the quadrupole and leave it 
by refraction preferrably near the points of highest cur- 
vature, (ii) The intensity loss per round trip should be 
neither too small nor too large. In the former case, too 
many individual orbits (with slightly different lengths) 
contribute to the far-field interference leading to an in- 
coherent response, whereas in the latter case the number 
of contributing beams is too small to produce a sharp 
interference pattern. This puts constraints on the refrac- 
tive index of the fiber. We have performed simulations 
with other fiber geometries and refraction indices which 
confirm the points (i) and (ii) above: using e.g. n = 1.5 
the filter effect disappears. 

Summary. Our ray-tracing model appears to describe 
the main features found in the experiment |3| . In particu- 
lar, the range of input and output angles, where far-field 
interference with filter characteristics can be observed, 



is rather small ("magic window"). The analysis of the 
length-difference spectrum allows for a clear distinction 
between our model of interfering rays from different or- 
bits and other scenarios involving a single orbit only. We 
note that the TE polarization used in Ref. J3| allows for a 
very efficient Brewster-angle input and output fiber cou- 
pling. 

Issues to be discussed in future work include (i) a com- 
parison to the case of TM polarization, (ii) the precise 
connection of the observed interference to chaotic cac- 
ity trajectories, and (iii) possible application of the filter 
effect for beam and/or cavity diagnostics. 
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